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EXPERIMENTAL AND ANALYTICAL DETERMINATION OF STABILITY 
PARAMETERS FOR A BALLOON TETHERED IN A WIND 

By L. Tracy Redd, Robert M. Bennett, and Samuel R. Bland 
Langley Research Center 

SUMMARY 


Experimental and analytical techniques for determining stability parameters for a 
balloon tethered in a steady wind are described. These techniques are applied to a par- 
ticular 7.64-meter-long balloon, and the results are presented. The stability parameters 
of interest appear as coefficients in linearized stability equations and are derived from 
the various forces and moments acting on the balloon. In several cases the results from 
the experimental and analytical techniques are compared and suggestions are given as to 
which techniques are the most practical means of determining values for the stability 
parameters. 


INTRODUCTION 


A tethered balloon is useful for a number of purposes, such as supporting antennas 
or providing an aerial platform. Such operations are often impaired by the occurrence of 
dynamic instabilities, especially during strong wind conditions. Although information 
relating to the stability of towed and tethered bodies, including balloons, has been pub- 
lished (see, e.g., refs. 1 to 4), a systematic procedure for the analysis of the stability of 
tethered balloons is needed. In an attempt to fill this need, the Langley Research Center 
undertook a general research study to develop improved techniques for predicting the 
stability of tethered balloons. 

Portions of the Langley study are given in references 5 to 8. Reference 5 describes 
the preliminary version of the towing technique that was developed and used to measure 
aerodynamic forces on inflated balloons. Reference 6 presents the derivation of the sta- 
bility equations, a comparison of analytical with experimental results, and a stability- 
parameter trend study. Reference 7 contains a listing and description of the computer . 
programs used for calculating and plotting the analytical results. A brief overview of 
the complete study is given in reference 8. 

The purpose of this paper is to describe methods of obtaining values for the sta- 
bility parameters required for the analysis of the stability of a tethered balloon. These 



methods are illustrated by their application to a particular 7. 64 -meter balloon used in 
the Langley stability studies. When possible, the analytically determined values of the 
parameters are compared with their measured values. 

SYMBOLS 


perpendicular distance from balloon axis of symmetry to a support 
cable as used in the roll-inertia test, m (see fig. 6) 


B 


distance along balloon center line from nose to structural center of 
mass, m (see fig. 2) 

buoyancy force, N 


b l’ b 2’ b 3 


'D 


'L,T 


distance along balloon center line from nose to center of buoyancy, m 
(see fig. 2) 

distances parallel to balloon center line from structural center of mass 
of balloon to front-cable attachment, rear -cable attachment, and 
center of buoyancy, respectively, m (see fig. 10) 


drag coefficient, 
lift coefficient, 


D 


pV 2 s/2 

L 


t 


pV 2 s/2 ] 

lift coefficient for the tail fins, 


pV 2 s/2 


K)t = 


9C l,t 

da 


rolling-moment coefficient, 




pV 2 Sc/2 


m 


pitching -moment coefficient, 




pV 2 Sc/2 


^m,T 

( Cm a)r 


Mv T’ 

pitching -moment coefficient for tail fins, — 


pV 2 Sc/2 


ac 


m,T 

da 


2 



yawing -moment coefficient, 


It 


M Z 

pV 2 Sc/2 


side-force coefficient, 


)V 2 s/2 


distance along balloon center line from nose to effective center of aero- 
dynamic apparent mass, m (see fig. 2) 


distances from balloon center of mass to front- and rear-cable attach- 
ment points, respectively, for yaw-inertia test, m (see fig. 8) 


balloon length, m 


aerodynamic drag force, N 

distance from balloon center of mass to elastic axis (e.a.), m (see 
fig. 8(b)) 


tensions in front and rear restraining cables used in aerostatic tests, N 
(see figs. 3 and 4) 

negative buoyancy force of inflation air for aerostatic tests, N 
force along the Y-axis, N 
acceleration of gravity, m/sec 2 

perpendicular distance from center line of balloon to load-band plane, m 
(see fig. 10) 

width of hull segment, m 

perpendicular distance from balloon center line to aerodynamic center 
of tail fin for planform A, m (see fig. 11) 

mass moment of inertia of compound -pendulum system about knife-edge 

O 

support, kg-m* 1 



w 


zz 


total moments of inertia about balloon center of mass, in the 

X",Y",Z" axis system (including aerodynamic apparent inertias), 
kg-m 2 (see fig. 2 for axis system) 


*xx a’*yy a’*zz a moments °f inertia of aerodynamic apparent mass about its own 


centroid, kg-m 


*xx,f 


roll inertia of trifilar -pendulum support frame, kg-m^ 


r 

l^x g>!yy g'^zz g moments of inertia of inflation gas about its own center of mass, kg-m^ 

Ixx s’*yy s’*zz s momen t s °f inertia of balloon structure about its own center of mass, 
kg-m 2 


A zz,f 

*zz,Q 


yaw inertia of the compound-pendulum support frame, kg-m^ (see fig. 7) 

combined yawing moment of inertia of the balloon structure, inflation 

o 

gas, and apparent mass about knife-edge support, kg-m (see 
figs. 7 and 13) 


K = m w^Q - B ( z Q - b ) + V( l Q " a ) + m f^f 


apparent-mass coefficient (see eq. (8)) 


k yy = 


T l + T 2 

h 



L 


l 2 


Tjlp + T n l 


2 R 


h 


aerodynamic lift force, N 

aerodynamic lift force due to tail-fin assembly, N 

length of tether cables used in bifilar -pendulum tests, m (see figs, 
and 8) 


length of support cables used in roll-inertia test, m (see fig. 6) 
distance from balloon nose to elastic axis (e.a.), m (see fig. 8) 
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distance from elastic axis (e.a.) to front tether cable for yaw-inertia 
test, .m (see fig. 8) 

distance from balloon nose to front tether cable as used in the yaw- 
inertia test, m (see fig. 8) 

distance from knife-edge support to center of gravity of support frame, 
m (see fig. 7) 

distance from knife-edge support to nose of balloon, m (see fig. 7) 

distance from elastic axis (e.a.) to rear tether cable for yaw-inertia 
test, m (see fig. 8) 

length of support frame for compound-pendulum test, m (see fig. 7) 

distance parallel to balloon center line from reference point to aero- 
dynamic center of tail-fin assembly, m (see fig. 11) 

distance between front and rear cables for yaw -inertia test, m (see 
fig. 8) 

distance along balloon center line from nose to center of mass of 
balloon (including aerodynamic apparent mass), m (see fig. 2) 

mean molecular weight of a gas, kg/mole 


aerodynamic rolling, pitching, and yawing moments about the X-, Y-, 
and Z-axes, respectively, N-m 

aerodynamic pitching moment due to tail-fin assembly, N-m 

mass of support frame used in compound-pendulum test, kg (see fig. 7) 

mass of inflation gas, kg 


W s 

mass of balloon structure without bridle or payload, — , kg 

g 



m w mass of ballast added to balloon nose for compound-pendulum test, kg 

(see fig. 7) 

in „,rn r o aerodynamic apparent mass associated with balloon acceleration along 

X,a. y,a Z,a 

X"-, Y"-, and Z"-axes, respectively, kg (see fig. 2 for axis system) 


m + ,m .,m , total mass of balloon for accelerations along X"-, Y"-, and Z"-axes, 

Xj L y y t Z y L 

respectively, for = 0, kg 
P oscillation period, sec/cycle 

p absolute pressure, N/m^ 


P,Q,r 


R* 

R 

S 

T 

T 1’ T 2’ T 3’ T 4 

fp ^ 

V 

V b 

W s 

X,Y,Z 
X”, Y", Z” 


rolling, pitching, and yawing rates about X-, Y-, and Z-axes, respec- 
tively, rad/sec 

universal gas constant, J /K-mole 

D * 

engineering gas constant, — , J /K-kg 

M 

reference area of balloon, V b 2// **, m 2 

sum of tensions in trifilar -pendulum cables, N (see eq. (11) and fig. 6) 

tensions in front and rear tether cables, N (see figs. 5, 8, and 10) 

absolute temperature, kelvins 

steady wind velocity, m/sec 

3 

volume of balloon hull, m 

structural weight of balloon without bridle or payload, N 

stability axis system; for steady-state conditions, X-axis is parallel to 
steady horizontal wind vector (see fig. 10 for axis system) 

body-fixed axis system; X"-axis is along balloon axis of symmetry (see 
fig. 2 for axis system) > 
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x,y,z 

x",y",z" 

X i 

x n 


y l’ y 2 


y 3 


displacements along X-, Y-, and Z-axes, m 

displacements along X"-, Y"-, and Z"-axes, m 

distance from balloon nose to center of mass of a hull segment, m 

distance from elastic axis (e.a.) to node line, m (see fig. 8) 

lateral displacements of upper ends of front and rear cables, respec- 
tively, for yaw-inertia test, m (see fig. 8(b)) 

lateral displacement of elastic axis (e.a.) in yaw-inertia test, m (see 
fig. 8(b)) 


a 


a t 

a lb 

“p 

0 




p 

p g 

<p 




<x> 


perturbation angle of attack of X-axis, z/V + 9, rad 
trim angle of attack, rad or deg 

angle of attack of tow-truck load bar, rad (see fig. 10) 

effective angle of attack of tail -fin assembly due to rolling rate p, rad 

angle of sideslip; rad 

angles of pitch, roll, and yaw, respectively, rad 
atmospheric density, kg/m^ 

density of inflation gas at standard atmospheric conditions, kg/m 3 

angular displacement of balloon tether lines in XZ -plane for bifilar- 
pendulum test, rad (see fig. 5) 

angular displacements of balloon tow lines in XZ-plane, rad (see 
fig. 10) 

frequency, rad/sec 
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Subscripts: 


air 


inflation gas is air 


g 

He 

t 

vac 


ambient air outside the balloon hull 


balloon inflation gas 


inflation gas is helium 


balloon at trim condition 


vacuum conditions 


a 

a 

P 

P 

P 

q 


Subscripts used with the aerodynamic coefficients denote derivatives as follows: 


with respect to a 


with respect to 
with respect to /3 
with respect to 


a c 
2V 


§c_ 

2V 


with respect to 2£. 
F 2V 

with respect to 

F 2V 


rc 

2V 


with respect to 

Dots over a symbol indicate derivatives with respect to time. 


GENERAL COMMENTS 

In the present study four types of parameters were required for the analysis of the 
stability of a tethered balloon: 

(1) Aerostatic properties (e.g., weight and buoyancy) 

(2) Masses and moments of inertia (including apparent mass) 
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(3) Aerodynamic coefficients and derivatives 

(4) Tether -cable properties 

These parameters appear as coefficients in the linearized stability equations of refer- 
ence 6 and are typical of those needed in other studies of the stability of towed or tethered 
bodies (refs. 1 to 4 and 9 and 10, e.g.). The values of these parameters, with the excep- 
tion of the tether-cable properties, were determined by techniques described in the sub- 
sequent sections of this paper. The tether-cable properties are not discussed in this 
paper since some of these properties were measured directly and the rest were deter- 
mined analytically by methods of reference 6. 

To obtain the values of the required balloon stability parameters, it was necessary 
to develop both analytical and experimental methods . The aerostatic properties were 
determined experimentally. The masses and moments of inertia were determined both 
analytically and experimentally, and the results were compared. The static aerodynamic 
lift, drag, and pitching -moment coefficients were determined experimentally from tow 
tests, and the remainder of the aerodynamic stability parameters were derived 
analytically. 

The 7.64-meter-long balloon used in the Langley stability study is shown in fig- 
ure 1, and its geometrical properties are given in table I and figure 2. 

The balloon geometry shown in table I was obtained by combining the forward por- 
tion of a C -class airship configuration (ref. 11) of fineness ratio 2 with a conical aft sec- 
tion to provide an overall fineness ratio of 3. The conical aft section joined the C-class 
portion at the point where the slopes were tangent. The nose of the balloon was also 
slightly more spherical than that of the regular C-class airship configuration. 

The essential components of the balloon were a hull (i.e., gas bag) constructed of 
nylon and saran, and a tail-fin assembly. A load band was attached to the lower portion 
of the hull to permit variations in the bridle attachment points and, hence, make it pos- 
sible to change the balloon trim angle of attack. The tail assembly consisted of four fins, 
each with a balsa spar framework covered with aluminized Mylar sheeting. The fins 
were connected together by rings at their front and rear and by a number of small cables 
that could be adjusted to aline the fins with the balloon center line. The tail assembly 
was connected to attachment, patches on the hull at four points for each fin. The balloon 
was designed to operate with a constant pressure differential; therefore, a relief valve 
was provided to avoid overinflation. 
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DESCRIPTION OF ANALYTICAL AND EXPERIMENTAL TECHNIQUES 


The balloon used in the present study was axially symmetric (fig. 1) except for the 
load band. Therefore, the effect of the load band was neglected and the masses, moments 
of inertia, and aerodynamic forces associated with motions in the z"-direction (fig. 2) 
were assumed to be the same as those for motions in the y"-direction. The values of 
the stability parameters thus obtained were for the axially symmetric balloon and did not 
include the effect of the bridle or payload. The effects of the bridle and payload (e.g., 
test instruments) were later added to the results of the present report to obtain the 
properties of the fully instrumented balloon as used for the study of reference 6. 

AEROSTATIC PROPERTIES 


The required aerostatic quantities a, b, B, nig, V b , and W s were determined 
experimentally from several tests. The basic structural weight W s was found by 
weighing the deflated balloon. The remaining parameters were determined from two dif- 
ferent tests with the inflated balloon constrained by two cables. In the first test (fig. 3) 
the balloon was inflated with helium; in the second test (fig. 4), with air. The first test 
was used to determine the buoyancy B from the measured tensions F 1 and F 2 by 
the equation 

B = Fj + F 2 + W g (1) 

The volume and mass of the inflation gas mg for the balloon were then computed 
by use of the perfect gas law (ref. 12) to give 


v b = 


B 


L _ P£ MjPe 

\ p e Rg 


R e T* 


( 2 ) 


m g = Pg V b (3 > 

where e refers to the external ambient air and g refers to the inflation gas. The 
volume can also be calculated from the balloon geometry. 

The locations of the structural center of mass a and of the center of buoyancy b 
are obtained by summing the moments about the nose for the two tests (figs. 3 and 4) to 
give 


aW g - bB = -cF 2 

r 


( 4 ) 
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and 


aW s + bFg = cF 4 


(5) 


The negative buoyancy F B of the inflation air in the second test was too small to mea- 
sure with the test instrumentation; therefore, it was computed from the perfect gas law 
as follows: 


F gV b(P e - Pe) 
B R e T* 


Solution of equations (4) and (5) yields 


c(bF 4 - F b F 2 ) 
W S (B + F b ) 


b = 


5 ( F 2 + F 4 ) 

B + F b 


( 6 ) 

(7) 


MASSES AND MOMENTS OF INERTIA 

In the present study both analytical and experimental methods were used to obtain 
the mass and inertial properties of the balloon. Some of the experimental methods were 
considered to be impractical for use with balloons larger than the 7. 64 -meter balloon. 
Therefore, analytical methods which could be used with larger balloons without difficulty 
were also developed. The validity of the analytical methods was evaluated by comparing 
the analytical results with those measured experimentally. 

Analytical Estimates 

The analytical techniques were applied by considering the balloon system to be 
composed of three parts: the inflation gas, the balloon structure, and the apparent air 
mass associated with the balloon accelerations. The mass and inertia of each of these 
parts were determined separately and then summed in an appropriate manner to obtain 
the mass and inertia of the entire balloon. 

Inflation gas . - The mass and inertial properties of the inflation gas were calculated 
by assuming the inflated balloon hull to be divided into 35 segments. These segments 
were formed by analytically slicing the hull normal to its axis of symmetry, as shown in 
sketch 1. The gas in each segment was assumed to be rigid and to form a frustum of a 
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right circular solid cone. The volume, mass, centroid location, and inertias about all 
three axes at the segment center of mass were calculated for each gas segment. The 
results were transferred from the segment centers of mass to the center of mass of the 
inflation gas and numerically integrated to obtain the total mass and inertial properties 
of the gas. 

Balloon structure .- The theoretical mass and inertial properties of the balloon 
structure were obtained by assuming the balloon to consist of two components: a hull 
and a tail-fin assembly. The properties of the total structure were obtained by appro- 
priately combining the properties of these two components. 

The hull mass and inertial properties were obtained by assuming the hull skin to be 
a rigid shell with a shape the same as that of the fully inflated balloon. The shell mass 
was assumed to be uniformly distributed (except for the seam tapes) and equal to the 
mass of the fabric of the balloon hull. The shell was divided into 35 segmental rings in 
the shape of truncated hollow cones, of the same widths hj as those used in calculating 
the mass of the inflation gas. The seam tape elements attached to each segment were 
assumed to be point masses appropriately distributed among the segments. The mass, 
inertias, and centroid location of each segment were calculated and the results were 
transferred to the shell center of mass. These results were then summed to give the 
mass and inertial properties for the total balloon hull. 

The mass and inertial properties of the tail-fin assembly were calculated by trans- 
ferring the mass and inertias of each tail fin and fin ring (fig. 1) to the assembly center 
of mass and summing the results. Similarly, the mass and inertias of each tail fin were 
determined by transferring the mass and inertias of each balsa spar and Mylar sheet to 
the fin center of mass and summing the results. 
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Apparent air mass . - When a body immersed in a fluid is accelerating, a certain 
amount of the surrounding fluid moves with the body. That portion of the fluid set in 
motion is referred to as the aerodynamic apparent mass, or simply apparent mass. 
Because the mass of the fluid (i.e., air) involved in the present study is of the same order 
of magnitude as the combined masses of the balloon structure and lifting gas, the apparent 
mass must be considered as part of the mass of the balloon system. 

The apparent mass of air associated with an elongated body accelerating parallel 
to its axis of symmetry is different from that for motions perpendicular to the axis of 
symmetry (ref. 13). Although the magnitudes of these masses are different, it is assumed 
in the present study that their respective mass centroids are all located on the axis of 
symmetry. Also, because of axial symmetry, it is assumed that the centroids are located 
at the same distance c (fig. 2) aft of the balloon nose for motions in both the y”- and 
z”-directions. 

The value of the apparent mass for acceleration along the X"-axis m x a was cal- 
culated by methods given in reference 13. In using these methods it was assumed that 
the aft end of the balloon hull came to a point and that the apparent mass associated with 
the tail fins for motions along the X"-axis was negligible. Furthermore, m x a was 
assumed to be proportional to V b and p; thus, 

m x,a = k lP V b ( 8 ) 

where kj is the apparent-mass coefficient which varies according to the width-length 
ratio of the hull (kj = 0.220 for the present analysis). 

The terms for the apparent masses and inertias, other than m x a , were calculated 
from procedures given in reference 14. In using these procedures, the balloon was 
divided into two sections, as shown in sketch 2, where section CD, or balloon "fuselage”' 
section, includes all the hull forward of the roots of the tail -fin leading edges; and 
section (2), or balloon "tail fin" section, includes the tail fins and the portion of the hull 
aft of the root of the tail-fin leading edges. 
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The equations for determining the apparent masses and inertias of the fuselage and 
tail-fin sections are given in detail in reference 14 and, therefore, are not presented in 
this report. In using these equations, it should be remembered that since the balloon is 
axially symmetric, m z a and I zz a are equal to m y a and Iy y a , respectively. The 
apparent roll inertia of the fuselage section is also assumed to be negligible when com- 
pared with that of the tail section. 


Experimental Measurements 

Four different types of test were required to measure the mass and inertial prop- 
erties. The first test was used to obtain thp apparent mass for motion along the X”-axis 
(fig. 2), the second was needed to measure the roll inertia, and the remaining two were 
used to determine the yaw inertia and mass associated with motions along the Y"-axis. 
Because the balloon was considered to be axially symmetric, the pitch inertia and mass 
associated with vertical motions along the Z"-axis were assumed to be equal, respec- 
tively, to the yaw inertia and mass associated with lateral motions. 

Longitudinal apparent mass m xa .- In order to determine experimentally the value 
of m x a , the balloon can be inflated with helium and tethered by two parallel cables to 
form an inverted bifilar-pendulum system, as shown in figure 5. If the longitudinal angle 
of oscillation <p of the parallel cables is small and the damping is neglected, the forces 
can be summed in the x-direction to give the equation of motion for the pendulum system: 


m x,t* + ( T l + T 2)i7 =0 
Since the period of oscillation P is 

P = 2tt 

and 


l T 1 + T 2 


m x,t = m s 


m x,a + 


m. 


then 


m 


x,a 


T 1 + T 2 p2 

4 ,% 



+ m g) 


(9) 


Since Tj, T 2 , and P are measured during the oscillation tests, and m g , and mg 
can be measured directly or by methods described in the section "Aerostatic Properties," 
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then m x a can be obtained directly from equation (9). (Note, eq. (9) can be ill condi- 
tioned if m x a « m g + 

Roll inertia.- In order to determine the roll inertia of a balloon, it can be suspended 
in a vacuum test chamber by three flexible cables of equal lengths and distances from the 
axis of symmetry to form a trifilar pendulum, as shown in figure 6. Assuming that the 
damping is negligible and the rolling oscillations are small, the equation of motion about 
the pendulum center line is 

The roll inertia 1^ can be expressed in terms of the period of oscillation P and also 
in terms of its component parts as follows: 

= *xx,a + *xx,f + *xx,g + *xx,s ( 10 ) 

The expression for the sum of the cable tensions T is found from the total weight and 
buoyancy of the pendulum system and is as follows: 


T = m g g + mjg - B 


( 11 ) 


If the oscillation tests are performed in a vacuum, then 1^ a and 1^ are 
negligible and the roll inertia of the balloon structure 1^ can be determined directly 
from equation (10), provided the roll inertia of the support frame is also known. Further- 
more, by subtracting the value of the roll inertia 1^ yac measured in a vacuum from 
the roll inertia 1^ obtained at atmospheric pressure, the sum of the roll inertias of 
the apparent mass and the inflation gas can be obtained from equation (10) as follows: 


Al xx *xx " *xx,vac 2, [^ P " ^ TP ^vacj *xx,a + *xx,g 

47T l 


(12) 


The combined values of 1^ a and 1^ g in equation (12) can be separated by 
testing the balloon at the same atmospheric pressure but with different inflation gases, 
for example, air and helium. An additional relationship derived from the perfect gas law 
(ref. 12) must also be used with equation (12). This relationship is 


(*xx,g)air ( I xx,g:) ; 


'He 


i RT *)airjp e=Const ,p 


( 13 ) 


where tile subscripts air and He designate air and helium inflation gases. 
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Yaw inertia and other mass properties .- The yawing moments of inertia 
l zz S’ l zz g’ and *zz ,a and the mass properties m y>t , m y>a , c, and Z x can be 
determined by the two pendulum tests shown in figures 7 and 8. For the first test (fig. 7) 
the balloon is suspended from a knife-edge support so that the balloon swings as a com- 
pound pendulum. For the second test (fig. 8) the balloon is inflated with helium and 
tethered by two cables to form a bifilar pendulum which has two degrees of freedom - 
yaw i// and side translation y. 

The compound-pendulum test (fig. 7) is a common experimental method for deter- 
mining the masses and moments of inertia of a body. (See, ref. 15, e.g.) In applying 
this method to the present case, special care must be taken to insure that the balloon 
retains its fully inflated shape during the tests. A ballast weight is also needed in some 
instances to provide sufficient restoring force to overcome the buoyancy of the inflation 
gas. 

The equation of motion for the compound-pendulum system used in the present 
study (fig. 7) is obtained by summing moments of inertia about the knife-edge support. 

If the damping is assumed negligible and the yaw angles i p are small, this equation is 

Iqi£+Kj//=0 (14) 


where 


K = m w^Q - B Pq - b ) + m s^Q " a ) + m f#f 
Thus the expression for Iq in terms of its period of oscillation P is 




KP^ 

4tt 2 


(15) 


The inertia Iq of the pendulum system can be expressed in terms of the masses 
and inertias of its component parts by use of the parallel-axis theorem to give 

>Q = hz,s + m s(^Q - a ) 2 + \z,t + m f ! f 2 + \z,t + “gpQ - b ) 2 

+ hz.n * m y,a(*Q * 0 f + m w*Q 2 < 16) 

where the moment of inertia of the ballast about its own center of gravity is assumed to 
be negligible. 

If the balloon is tested in a vacuum, then the inflation-gas and apparent-mass terms 
in equation (16) become negligibly small. Thus, for these conditions the inertia of the 
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balloon structure Iyy s can be determined from equation (15) and the remaining expres- 
sions in equation (16), provided the mass and inertia of the support frame and ballast 
weight are known. The mass and inertia of the support frame can, of course, be deter- 
mined by oscillating the frame without the balloon. By subtracting the inertia Iq vac 
obtained in a vacuum from the inertia Iq measured at atmospheric conditions, the 
inertias of the apparent mass and the inflation gas are obtained as follows: 


- ’Q - *Q,vac - [ Kp2 - ( Kp2 )vac]p 

where 

A] Q = W + m gfo ‘ b ) 2 + W + m y,a?Q - C ) 2 


(17) 

(18) 



and w is the oscillation frequency of the system. 

Equations (18) to (23) provide a set of six equations with seven unknowns, namely, 
I zz , l zz a , l zz g, m yt , m y a , l x , and c. All other terms in these equations can be 
obtained from direct measurements or by tests described in the section entitled "Aero- 
static Properties" and in the appendix. 

Since the preceding equations contain one too many unknowns, it is suggested that 
the analytical value of Iyy )g be used because this term is the most straightforward to 
calculate. Thus, equations (18) to (23) can be solved simultaneously to find the six 
remaining unknowns. 
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AERODYNAMIC COEFFICIENTS AND DERIVATIVES 


Lift, Drag, and Pitching -Moment Coefficients 

The aerodynamic forces on an inflated balloon are generally obtained from wind- 
tunnel tests of a rigid model (refs. 16 and 17, e.g.). Such test methods do not account for 
deformations of the flexible balloon structure, which may be significant at certain load 
conditions. For this reason, a tow -test technique was developed to measure the static 
aerodynamic forces and moments acting on a flexible inflated balloon at various angles of 
attack and wind speeds. 

The basic difference between the tow -test technique and conventional wind-tunnel 
methods is that the forces and moments on the balloon are measured from the tensions in 
the tow cables rather than on conventional force balances. Thus, the towed balloon does 
not have to be rigidly mounted and can have structural deformations and some freedom of 
motion. 

Description of test technique .- The tow-test technique used in this study is illus- 
trated in figures 9 and 10. In figure 9, the balloon is shown attached to a panel truck by 
two cables and a load bar. The panel truck provides space for instrumentation and acts 
as a towing vehicle. The load bar is attached to the truck by means of a pivot so that its 
angle of attack can be varied. The load bar is also equipped with load cells to measure 
the cable tensions and with potentiometers to measure the cable angles with respect to 
the load bar. The speed of the truck relative to the air and the sideslip angle are mea- 
sured by an anemometer (i.e., wind sensor) mounted at the front of the truck. Outputs 
from these instruments are recorded on strip charts. Motions of the balloon are also 
recorded by' two cameras mounted on the truck. 

Determination of aerodynamic coefficients .- When the system (fig. 10) is in equilib- 
rium, the aerodynamic forces and pitching moment L, D, and My on the balloon are 
balanced by the balloon weight W g , buoyancy B, and cable tensions T3 and T4. For 
this condition, the expressions for the lift, drag, and pitching moment about the refer- 
ence point (figs. 2 and 10) are 

L = W s - B + T 3 cos <P 3 + T 4 cos cp 4 (24) 

D = T3 sin cp 3 + sin cp 4 (25) 

My = -b3B cos + TgJbj cos^g - - h sin(<£>3 - 

+ T 4 |b 2 cos (<p 4 - a t ) - h sin(<p 4 - a? t )j (26) 


18 



The quantities W s , bj, b 2 , b3, and h in these equations are determined either from 
direct measurements or from the aerostatic tests described previously. The quan- 
tities a t , T3, T4, cp 2 , and can also be measured directly in some instances, 
but for the present configuration it is more convenient to calculate these quantities from 
the measured values of Tj, T 2 , <P\, and cp 2 (fig. 10 ) at the lower ends of the 

cables. Thus for this configuration, the additional steps which are explained in the next 
few paragraphs are required to obtain the lift, drag, and pitching moment of the balloon. 

At zero wind speed, the balloon, load bar, and tow cables form a parallelogram, so 
that the balloon is at the same angle of attack as the load bar and the cables are vertical. 
With the truck moving, the drag of the balloon causes it to ride behind its former posi- 
tion, so that the cables are no longer vertical. In this position, the drag on the cables 
and their weight act to bow the cables. In addition, if the tensions in the two cables are 
unequal, they will bow by different amounts, so that the balloon is no longer parallel to 
the load bar, as shown in sketch 3 . 



To relate the tension and angles at the upper ends of the cables to those measured 
at the lower ends, the analysis described in appendix B of reference 6 is used. This 
analysis yields the position, tension, and angle of the upper end of a cable in terms of the 
cable length, weight, diameter, drag coefficient, wind speed, and the tension and angle at 
the lower end of the cable. For the present application, this is done for each cable inde- 
pendently; at first the constraint that the upper ends of the two cables are actually a 
known, fixed distance apart is ignored. This constraint is then taken into account by 
adjusting the measured angles at the lower ends of the cables by a small amount and 
repeating the cycle. This process would not have to be repeated if the cable analysis 
(ref. 6) gave exact results and if the cable angles cpi and <p 2 and tensions T4 
and T 2 could be measured precisely. Since this is not the case in the present study, 
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it is desirable to repeat the process in an iterative fashion until the predicted distance 
between the upper ends of the cables agrees with the known value. 

Equations for the lift, drag, and pitching -moment coefficients can be obtained by 
plotting Cl, Cq, and C m as a function of a ^ and fitting least-squares curves 
through the resulting data points. These equations can be differentiated with respect 

to a + to obtain expressions for C T , Cn , and Cm . 

1 a a a 

Calculation of Longitudinal Dynamic Stability Derivatives 

The stability derivatives C T , C™ , C T . , and C™ . are needed in the sta- 

Uq niq m a 

bility analysis (ref. 6) in addition to the longitudinal static coefficients that were mea- 
sured. In the present study these dynamic derivatives are determined analytically for 
use in the stability analysis. 

Hull contributions .- The contributions of the balloon hull to the dynamic derivatives 
could not be analytically determined. The difficulty was that the static aerodynamic data 
measured for the hull with and without tail fins were insufficient to use nonslender body 
theory (e.g., ref. 18) for calculating these contributions. The values of CLq and C m q 

were obtained by using slender-body theory (e.g., refs. 18 and 19) and were found to be 
small compared with the tail-fin contributions. On the basis of these slender-body 
results and in the absence of further data, the hull contributions for the present configu- ' 
ration were neglected, even though the hull is not a slender body in the usual sense. 

Tail -fin contributions .- Both static and rate stability derivatives were calculated 
for an isolated, planar tail fin by applying linear, oscillatory, lifting -surface theory to 
downwash deformation modes of camber, pitch, and plunge, and deriving the stability 
derivatives from the resulting generalized aerodynamic forces in the manner discussed, 
for example, in reference 20. The numerical method used, which was written and pro- 
gramed by Robert N. Desmarais of Langley Research Center, is an unpublished develop- 
ment of the method described in reference 21. On each semispan, 36 collocation points 
were used in the calculations. 

One question that arises in the calculation of the tail-fin contributions is the treat- 
ment of the effective planform for the low-aspect-ratio surfaces attached to the relatively 
large, conical aft portion of the present balloon. Because of the uncertainties involved, 
three alternate ways of treating the effective planform were considered in the present 
study, as illustrated in figure 11. For planform A, the conical portion of the balloon is 
treated as a reflection plane and the free stream is considered to be parallel to the coni- 
cal portion; for planform B, the conical aft portion of the body is included rearward from 
the intersection of the leading edge of the fins with the body; and for planform C, the 
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leading edges of the fins are projected to the center line of the body. Of the three plan- 
forms evaluated, planform A is considered to be the most realistic. The reasons for 
this conclusion are given in the section entitled "Dynamic Stability Derivatives." 

Estimates of Lateral Stability Derivatives 

None of the lateral stability derivatives were measured; hence, all these deriva- 
tives had to be estimated by other means. As was pointed out earlier, the balloon con- 
figuration was assumed to be axially symmetric. Thus, for small values of a the 
following lateral derivatives are assumed to correspond to the longitudinal derivatives 
for = 0: 



_CL «L t =o 

Cn |3 = 

" Cm «L t =0 

c Yr 

“ CL aL t =:0 

Cn /3 = 

' Cm «L t =0 

C Y r = 

% 

o?t = 0 

c n r = 

Cm q 

ff t =° 


Note that the lateral derivatives are based on the length of the balloon c, not the span. 

Since the balloon is assumed axially symmetric, the lateral aerodynamic center is 
assumed to be on the center line. Thus, when a ^ = 0, the rolling moments due to side- 
slip /3 and yawing motions r are zero ^i.e., = 0 for = oj. It is 

also assumed that C, . = 0 for small values of a. 

If it is assumed that the lateral aerodynamic center of the balloon remains on the 
center line but displaces with angle of attack as the trim angle changes from zero, 
then Cj „ becomes 



sin ofj. = -C 


m a 


sin 
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The derivative C i r for a ^ different from zero can be approximated by neglect- 
ing the aerodynamic forces on the hull. For this case the Y-force due to yawing velocity 
is assumed to act on the balloon center line at the position of the static aerodynamic cen- 
ter of the tail -fin assembly, a distance (fig. 11) aft of the reference point, and dis- 
placed below the X stability axis by an amount equal to sin Thus becomes 


Zrp 

C 7 = -Cv -4- sin 
l v ±r c 1 


The effective angle of attack <Xp of the tail -fin assembly due to rolling rate 
(ref. 22) 


l 

p is 


a . 


= f *T sin «t 


Thus, if the hull effects are neglected, then the side -force and yawing-moment coefficients 
due to rolling motions are 


C Y= ( C h a ) T a V 

C n = " C Y Y cos <*t 

where (Cl^ is the lift-curve slope calculated for the 
Hence, the derivatives of C Y and C n with respect to 



The damping-in-roll derivative Cb was calculated in a similar fashion to that 

r 

for the longitudinal derivatives . For planforms B and C the derivative was derived 

directly from rolling about the center line of the balloon at = 0. For planform A the 
rolling about the root chord was treated and the resulting rolling moment transferred to 
the center line of the balloon. 


planform A tail-fin assembly, 
p are 
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It is recognized that these estimates of the lateral derivatives are simplified and 
may be somewhat crude. For example, estimates of C n for airplane-type configura- 

ir 

tions involve yawing moments generated by the change in induced drag due to rolling 
velocity (ref. 22 ,), which have been omitted here. Further detailed treatment was con- 
sidered unwarranted, however, because of the absence of body-fin interaction data. 

TEST PROCEDURES AND DISCUSSION OF RESULTS 
AEROSTATIC PROPERTIES 

The measured aerostatic properties defined by equations (1) to (7) and illustrated 
in figures 3 and 4 are shown in table II. Each of the aerostatic tests was performed at 
least three times and the data were averaged to obtain the results shown in the table. 

The scatter in the data was less than 1 percent except for and B, which varied as 
much as 1.3 and 2.4 percent, respectively. Most of this scatter in the latter two param- 
eters is attributed to changes in the temperature of the ambient air. 

MASSES AND MOMENTS OF INERTIA 

The results of the mass and inertia tests are shown in table HI. This table gives 
the calculated and measured values of the masses and inertias of the balloon and its 
component parts. Since some of the experimental data are believed to be in error, the 
table also includes a third column, "Best estimate," which combines both experimental 
and calculated data to obtain the most reasonable values of the actual balloon masses and 
inertias. The values of the terms my^, m z> t, Iyy, and I zz in the "Best estimate" 
column are each obtained by combining several balloon properties, some analytically 
determined and some experimentally determined. Thus, the best-estimate values for 
these terms are different from either their analytical or experimental values. The 
particular combination of balloon properties used to obtain the best-estimate values of 
these terms will be presented later in this section. 

All the data shown in table HI are for the axially symmetric balloon inflated with 
helium and surrounded by air at standard atmospheric conditions. 

Longitudinal Apparent Mass 

The experimental and calculated values of the apparent mass associated with longi- 
tudinal motion m x a are given in table HI. The exact agreement between these two 
values is believed to be a coincidence since the measured value could only be determined 
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to an accuracy of about 3 percent and the apparent mass of the tail fins was neglected for 
the calculated value. 

The quantities T^, T 2 , and P in equation (9) were measured with load links and 
potentiometers located at the base of the tether cables (fig. 5). The tests were repeated 
three times and the results were averaged to obtain the value of m x a shown in table HI. 

Roll Inertia 

The trifilar -pendulum test (fig. 6) used to obtain the roll inertia of the 7.64-meter 
balloon was conducted in the 20-meter-diameter vacuum sphere at the Langley Research 
Center at air densities ranging from 0.012 to 1.29 kg/m^. The periods of the roll oscil- 
lations were measured with a rate gyro and recorded on a strip recorder. 

The combined roll inertias of the balloon structure, inflation gas, and apparent 
mass are plotted as a function of the vacuum -chamber density in figure 12. This figure 
shows that for the air -inflated balloon, the inertia of the apparent mass and inflation gas 
is a linear function of the test-chamber density. Test points at chamber densities other 
than at vacuum and atmospheric conditions were obtained for this case so that a best-fit 
straight line could be drawn through the points to minimize the effect of scatter in the 
data. This procedure gave more accurate values for 1^ at vacuum and atmospheric 
conditions . 

The helium -inflated balloon was tested at chamber densities of 0.012 and 1.26 kg/m3. 
For the higher density, the balloon was tethered from the vacuum -chamber floor instead 
of the ceiling (fig. 6) because the buoyancy force provided the restoring moment. For the 
lower density, the balloon was suspended from the ceiling and the structural and frame 
weights provided the restoring moment. There were no tests for chamber densities 
between these two points because the combination of buoyancy, structural weight, and 
frame weight did not provide sufficient restoring moments for good oscillatory motions. 
Although data were available for only the two chamber densities, it was assumed that the 
inertia of the apparent mass and inflation gas for the helium -inflated balloon was a linear 
function of chamber density. Thus, a best-fit straight line was also drawn through these 
data to minimize scatter. 

The values of the roll inertias obtained from the vacuum -sphere tests and from the 
analysis are given in table III. This table shows that the analytical and measured values 
of Ijqj s are in good agreement, but that the experimental values of Ixx,g anc * ^xx,a 
are considerably smaller than those calculated. The calculated value of g is 
believed to be quite accurate when all the inflation gas rotates with the hull. For the 
present case, however, it appears that only the portion of the gas near the hull skin 
rotates with the balloon. Therefore, it is believed that the measured value of 1^ g is 
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the more realistic. The measured value of is also believed to be more accurate 

than the calculated value since the roll-inertia test is a straightforward test and some 
empirical estimates (ref. 14) were required to obtain the calculated value. 

Yaw Inertias and Other Mass Properties 

The desired yaw inertias and mass properties I zz , I zz s , I zz a , m y a , c, 
and l x for the 7. 64 -meter balloon were very difficult to determine experimentally 
because these properties were found to be functions of small differences between large 
values. Thus, in order to improve the prospects of obtaining accurate experimental 
values with the present approach, four different tests were performed with the compound- 
pendulum system (fig. 7). Different support lengths lg were used for three of the tests 
and the fourth test used a different inflation gas. The results from each of these tests 
were independently combined with the results from the bifilar test (fig. 8) to give four 
separate values for each of the desired mass properties. These values for each property 
were then averaged to obtain the final experimental results listed in table El. 

To increase the accuracy of the compound-pendulum test results further, enough 
data points were obtained to plot the yawing moment of inertia as a function of chamber 
density for each of the four tests. Since the yaw inertia was assumed to be a linear func- 
tion of density (ref. 15), a best-fit straight line was drawn through the points for each test 
to minimize the effect of scatter in the data. Other approaches similar to those of refer- 
ence 15 were attempted in order to obtain the desired mass properties, but the results 
were much worse than those obtained with the present method. 

The compound-pendulum tests were conducted in the 20-meter-diameter vacuum 
sphere at the Langley Research Center. The test-chamber air densities ranged from 
0.0016 to 1.26 kg/m^. The three different support -frame lengths lg used were 3.17, 
6.24, and 9.24 meters, and the inflation gases were air and helium. A ballast of 4.48 kg 
was attached to the nose of the balloon in order to maintain a sufficient restoring force 
when the balloon was filled with helium. The periods of oscillation for the compound- 
pendulum tests were measured with an accelerometer and recorded on a strip-chart 
recorder. 

The results of the compound -pendulum tests are shown in figure 13. The data pre- 
sented are for the combined inertias of the balloon structure, inflation gas, and apparent 
air mass about the knife-edge support (fig, 7) as a function of vacuum -chamber density. 

As can be seen from figure 13, the inertias of the apparent mass and inflation gases are 
proportional to the chamber density, as expected. Each data point shown represents the 
average of three or more compound-pendulum test results. The scatter at each of these 
data points is less than 1 percent; hence, individual test results are not shown in the 
figure. 
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The helium inflation gas was used only for the longest support frame length ig 
(fig. 13) in order to economize on time during the compound-pendulum tests. Further- 
more, the helium -inflated balloon configuration could not be tested in vacuum -chamber 
densities greater than 0.608 kg/m® because of the large buoyancy force at the higher 
densities. Additional ballast could have been used, but this would have deformed the 
nose appreciably. 

Some of the experimental values shown in table in are considerably different from 
the corresponding calculated values. The cause of these differences is attributed to 
scatter in the experimental data. Although the data for q (fig. 13) have less than 
1 percent scatter, this small amount of error combined with the errors from the bifilar- 
pendulum tests caused large errors (greater than 80 percent in some cases) in the exper- 
imentally determined values of I zz , Izz,s> Izz,a> m y,t> c > an< * lx- The major prob- 
lem was that the support -frame lengths £g were too long and thus caused the terms in 
equations (16), (18), and (22), such as m s (ZQ - a) 2 , m g (ZQ - b) 2 , Iq, and AIq, to be 
much larger than I^g or I zz>a . Hence, a small error in a, b, Iq, m s , m g , 

AI Q’ or Iq resulted in large errors in I^, l zz s , and I zz a . As a consequence, the 
values of my a , c, and l x were also somewhat questionable. 

Because of the uncertainty in the experimental results, the analytical values of 
I zz s’ *zz,a> m y,a> c » were believed to be more accurate. For the same rea- 

son, the best-estimate values for l zz and my ^ were not obtained from the experi- 
mental data but from equations (21) and (22) by combining the appropriate quantities from 
the "Best estimate" column of table in. 


AERODYNAMIC COEFFICIENTS AND DERIVATIVES 


Lift, Drag, and Pitching-Moment Coefficients 


The tow tests to determine the aerodynamic lift, drag, and pitching-moment coeffi- 
cients of the 7.64-meter balloon (figs. 9 and 10) were performed on a runway taxi strip 
at the Langley Research Center. The balloon was towed during periods of calm air, at 
constant speeds of from 3 to 15 m/sec (corresponding to Reynolds numbers based on c 
of from about 1.5 x 10® to 7.5 x 10®) for an angle-of-attack range from 0° to 24°. 

A typical strip-chart record of the tow-test measurements is shown in figure 14. 
As is evident in this figure, some of the data showed appreciable oscillations; hence, it 
was necessary to fair through the data to obtain the steady -state values of the aerody- 
namic parameters. 


In order to solve equations (24) to (26) for the lift, drag, and pitching moment, it 
was necessary to obtain the terms oe^, cp^, 


Tg, and 


T^ from the measured 


r' 
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quantities <pj, Tj., and T 2 by the iterative procedure discussed previously. 

For the present tests, however, it was necessary to limit the procedure to one iteration 
for two reasons. First, the adjustment in the measured angles cp^ and <p 2 had to be 
kept within what was believed to be a reasonable experimental error (3°). Second, 
although the procedure converged for small angles of attack, it did not appear to do so at 
larger angles. 

There was appreciable scatter in the data for repeated tests at a given angle of 
attack and velocity. Part of this scatter, as mentioned in reference 5, was due to small 
variations of ambient wind, limited length of the test runway, and oscillations in the 
traces on the strip-chart records. A large portion of the scatter was also due to changes 
in buoyancy caused by variations in the ambient temperature of the atmosphere. These 
variations led to uncertainties in the exact magnitude of the buoyancy force during many 
of the tests. To help overcome this problem, the aerodynamic forces and pitching 
moment (L, D, and My) were plotted against dynamic pressure pV^/2 at each angle 
of attack, and least-squares straight lines were fitted to the resulting data points. The 
aerodynamic coefficients for each angle were then derived by taking the slopes of these 
straight lines. The use of this method eliminated the need for the absolute value of the 
buoyancy force for every test and, hence, reduced the scatter in the aerodynamic coeffi- 
cients from that given in reference 5. There were not enough test data to use this 
method for angles of attack above 16° or for tests in which the tail fins were removed; 
thus, no data are presented for these cases. 

The lift, drag, and pitching-moment coefficients C L , Cq, and C m for the 
7.64-meter balloon are shown in figure 15. The curves shown are least-squares fits to 
the data and are given by 

C L = 0.82(o! t - 0.023) - 5.023(a t - 0.023) 3 + 111.4(a t - 0.023) 5 

C D = 0.0487 + 186.2 (a t - 0.023) 6 

C m = -0.0106 - 0.1435ce t 

where ce^ is the angle of attack in radians. The particular degree of the polynomials 
chosen for each of these equations is a matter of judgment; however, the lift should be an 
odd function and the drag an even function of angle of attack, as given. In the present 
case, the curves are weighted so that the data points close to = 0 have the most 
influence on the shape of the curves. Note that the angle of zero lift (and minimum drag) 
is taken to be = 0.023 radian (1.3°). 
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The shapes of the lift, drag, and pitching -moment curves (fig. 15) are somewhat 
different than expected for a normal C-class balloon. (See, e.g., refs. 16 and 17.) For 
the C-class balloon, the lift curve is generally more linear and the slope is steeper at 
lower angles of attack. Also, the drag curve generally starts to increase at lower angles 
of attack. These differences in the shapes of the curves are believed to be due primarily 
to detailed differences in the shapes of the C-class and present balloon configurations, 
particularly in the tail -fin assembly. The C-class tail -fin assembly generally consists 
of three thick fins which are not in the form of a delta wing (refs. 16 and 17), whereas the 
present configuration has four thin delta-wing fins. The lift-curve slope for the present 
case does increase with increasing angle of attack, which is consistent with the lift data 
for delta wings of low aspect ratio (ref. 23). 

Some of the differences in the shapes of the lift, drag, and pitching-moment curves 
for the C-class balloon and the present balloon are also believed to be caused by defor- 
mations of the nose and tail fins of the present balloon. In an extreme case (e.g., 

V = 15 m/sec and = 23°), the horizontal tail fins deformed to the extent that they 
failed structurally. 

From the results of this investigation, it would appear that the present tow-cable 
technique is adequate for obtaining the lift, drag, and pitching moment for inflatable bal- 
loons. This technique is particularly useful when deformations of the balloon structure 
or shape are expected. If the balloon structural deformations are not important, then 
wind-tunnel tests with rigid models may have some advantage over the present towing 
technique. For rigid models, the data are not contaminated by buoyancy changes and the 
static lateral aerodynamic forces and moments can also be measured. 

Dynamic Stability Derivatives 

The results of the calculations of the dynamic stability derivatives for each plan- 
form shown in figure 11 are given in table IV. These results are based on reference 
values of c and S for the complete balloon and include only the tail -fin contributions 
to the derivatives. 

As previously noted, the calculated stability derivatives for planform A are believed 
to be the most realistic of those presented in table IV. There are three reasons for this 
conclusion: 

(1) The calculated value of (Cj^ for planform A is nearer the measured value 
of Clq, for the complete balloon. Note that the calculated value of (CL a ) T for each 
of the rigid tail -fin planforms presented in table IV is higher than the measured value for 
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the complete balloon (measured Ct, J = 0.82). This condition could possibly result 

\ a| «t=0 / 

from aeroelastic deformations of the balloon. 

(2) The calculated value of (c ma ) T for planform A is nearer the measured incre- 
ment of Cm a due to the tail fins. Although the experimental data measured with the 
tail fins on and off the hull are inadequate for a quantitative definition of tail -fin effec- 
tiveness, there were enough tail-off data to get a qualitative estimate of the tail-fin con- 
tributions to the measured value of C ma ,. 

(3) Photographs of tufts on the lower surfaces of the horizontal tail fins indicate 
that except near the tips, the flow is nearly parallel to the conical aft portion of the body, 
as assumed for planform A. 

The values of C m q estimated by calculations are quite close for all three plan- 
form treatments (table IV). Larger differences due to planform treatment are apparent 
in the a derivatives, which generally arise from downwash-lag effects. The principal 
source for downwash-lag effects for conventional airplanes is the lag of downwash from 
the wing impinging on the tail (e.g., ref. 22). For the balloon configuration treated herein, 
interactions between the body and tail fins might be a similar source of downwash-lag 
effects, but these interactions could not be determined for the present case. Thus, these 
interactions are omitted here and the calculated values of and ^ m Qr take into 

account only the effects of the fin on itself. The calculated value of C m ^ may be under- 
estimated as a result of the omission of body-fin interaction effects. The values of CLq 
and Q p for planform A are quite different from those for planforms B and C. These 
differences are due to differences in planform areas and locations of the aerodynamic 
center. 

It should be noted that the linearized lifting- surface theory used to calculate the 
aerodynamic derivatives treats only the conditions for small angles of attack and does 
not include such effects as leading-edge vortex formation on highly swept planforms at 
angle of attack (e.g., refs. 23 and 24). 

CONCLUDING REMARKS 


Analytical and experimental methods have been presented for determining the sta- 
bility parameters for an inflatable balloon. These methods were applied to a 7.64-meter- 
long balloon, and the results were presented. 

Three types of stability parameters were considered in the present study. These 
were aerostatic properties, masses and moments of inertia, and aerodynamic coefficients 
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and derivatives. (The tether-cable derivatives were not included in this study but are 
presented in NASA TN D-7272.) 

The masses and moments of inertia for the 7. 64 -meter balloon were determined 
from both analytical and experimental methods and the results of the two methods were 
compared. Most of the analytical results were consistent with those obtained experi- 
mentally. There were significant differences, however, between the measured and calcu- 
lated values of apparent-mass inertias. These differences were attributed in some cases 
to the use of inadequate empirical estimates in the analytical calculations and in other 
cases to scatter in the experimental data. 

The results of the present study indicate that the static aerodynamic lift, drag, and 
pitching-moment coefficients for inflated balloons can be satisfactorily measured with a 
double-cable towing technique. This procedure is useful when there is difficulty in simu- 
lating elastic deformations of the balloon in a wind tunnel. 

The values of the aerodynamic stability derivatives calculated in this report may 
be somewhat crude, but they are believed to be of sufficient accuracy to warrant use in a 
tethered -balloon stability analysis. 

Langley Research Center, 

National Aeronautics and Space Administration, 

Hampton, Va., April 25, 1973. 
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APPENDIX 


EQUATIONS OF MOTION FOR THE BIFILAR PENDULUM 


The equations of motion for the bifilar -pendulum system (fig. 8) are derived by 
summing forces in the y-direction and moments about the Z"-axis (i.e., yawing moments). 
If the motions in the y- and ^-directions are small and the damping is negligible, then 
the equations of motion about the balloon center of mass are 


I 

I 


F v = m. 


y,t y + 


: . T iyi . ^2 


h 


= 0 


T l c l 


M z = I zz^ /+ ~ — yi + 


T 2 c 2 


y 2 = 0 


(Al) 

(A2) 


The point labeled "e.a." in figure 8 is at the intersection of the balloon center line 
and the elastic axis. The elastic axis is defined as the locus of points at which a lateral 
force can be applied and the balloon will translate laterally but will not rotate. Hence, 
about the elastic axis the following relationship applies: 


T 2 Z R =Ti* f 


(A3) 


If the origin of the coordinate system is transferred to the point "e.a.,” then the 
following relationships apply (fig. 8(b)): 


y i = y 3 + 1 f^ 

y 2 = y 3 - 1 r * 

y = y 3 - ej* 

C 1 = e i + Zp 

c 2 = e l " l R 
e l = l x~ l e 
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APPENDIX - Continued 


and the equations of motion (Al) and (A2) become 


where 


m y,t^3 + k yy y 3 " m y,t(^x “ l e)^ = 0 

l zz^ + k ^ +k yy{ l x- l e)y3 =0 


k yy = 




T l + T 2 


T ^ F 2 + T 2^R 2 


(A4) 

(A5) 

(A6) 

(A7) 


The expression for the location of the vertical axis which has no side translation 
(i.e., node line) is as follows (fig. 8(b)): 

x n =^ 

n ijs 

Substituting this relationship into the equations of motion (A4) and (A5) gives 

K - •• 


™ n x e 

y.tl x, 




I zz^ + ^ +k yy( i x- Z e) x n ]^ =0 


(A8) 

(A9) 


Since the damping is assumed negligible, the solution to equations (A8) and (A9) is 
in the form of 


^ 0 e 


icut 


where is a constant. Substituting this relationship into equations (A8) and (A9) 
gives 


m M. _ lx " l * \ = 

y-n x n J w 2 


i -a _ i w = 
zz l‘x L eJ n 2 2 

c o u> 


(A10) 

(AH) 
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APPENDIX - Concluded 


The system analyzed has two degrees of freedom and therefore will have two pos- 
sible modes. Each mode will have its own unique frequency u> and corresponding 
node-line location x n . In most cases, however, the system can be excited so that one 
mode will dominate the motion. Thus for this condition, the system may be considered 
to exhibit only one mode of motion. 

The terms l e , ky y , and in equations (A10) and (All) can be calculated 

from equations (A3), (A6), (A7), and the relationships (see fig. 8(a)) 

l R + l F = Z TH 
l e = l FC + l F 

The frequency co and the location of the node line Xjj can be determined directly from 
the bifilar- pendulum test. Thus, the remaining unknowns in equations (A10) and (All) 
are l^z’ ^y,t* ^x - 
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Distance along balloon 


Station axis of symmetry Radius, m 

from nose, m 


0 0 0 

1 .05 .25 

2 .10 .38 

3 .20 .60 

4 .41 .88 

5 .61 1.06 

6 .81 1.15 


7 1.02 1.21 

8 1.22 1.24 

9 1.42 1.27 

10 1.63 1.27 

11 1.83 1.27 

12 2.03 1.26 

13 2.24 1.24 

14 2.44 1.22 

15 2.64 1.19 

16 2.84 1.16 

17 3.05 1.12 

*18 3.25 1.08 

19 4.88 .71 

20 7.64 .10 










TABLE E.- AEROSTATIC PROPERTIES OF THE 7.64-METER BALLOON 

Buoyancy, B, N 190 

Center of buoyancy, b, m 2.48 

Mass of inflation gas, nig, kg 3.22 

Reference length, c, m 7.64 

Structural center of mass, a, m 4.72 

Structural weight, W s , N 84.2 

Volume, Vfo, m^ 19.0 
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TABLE DI.- ANALYTICAL AND EXPERIMENTAL VALUES OF MASS 
AND INERTIA FOR THE 7.64-METER BALLOON 


Parameter Analytical Experimental Best estimate 


Structure 


a, m 4.51 4.72 4.72 

m s , kg 8.28 8.58 8.58 

Ijoj s , kg-m^ 6.64 6.81 6.81 

I yy,s = ^zz,s> kg-m^ 56.5 54 56.5 


Inflation gas 


b, m 2.65 2.58 2.58 

mg, kg 3.08 3.22 3.22 

Ixx g, kg-m^ 1.83 0.24 0.24 

Iyy,g ~ I Z z,g> kg-m^ 7.90 7.90 


Apparent mass 


c, m 3.56 2.79 3.56 

m x a ,kg 5.11 5.11 5.11 

m y,a =m za , R g 23.9 23.0 23.9 

r xx,a> kg-m 2 2.78 1.58 1.58 

*yy,a = ^zz,a» k g _m2 76.5 19 76.5 


Total balloon* 


l x , m 3.70 3.25 3.70 

m x t> kg 14.3 16.9 16.9 

m y>t = m z t , kg 35.3 34.8 35.7 

1^, kg- m 2 11.2 8.63 8.63 

Iyy = ^zz> kg-m 2 150 105 155 


*The magnitudes of the terms l x , m x m y m z 1^, I yy , and l zz are 
different in reference 6 because the mass of the bridle and attached test instruments is 
added to the mass of the hull structure. 
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TABLE IV.- CALCULATED TAIL-FIN CONTRIBUTIONS 


TO STABILITY DERIVATIVES* 


Derivative 

Tail 

-fin contributions for - 

Planform A 

Planform B 

Planform C 

( C ^) T 

0.952 

1.19 

1.21 


.089 

.241 

.342 

% 

.685 

.847 

.855 

( Cm a) T 

-.247 

-.236 

-.182 

Cm a 

-.026 

-.061 

-.066 

c m 

m q 

-.189 

-.201 

-.193 

c , 

£ p 

-.0237 

-.0156 

-.0156, 


$ 

Based on S and c for complete balloon; moments are about the 
reference point. 
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Figure 1.- The 7.64-meter balloon. 
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Figure 2.- Geometrical properties of the 7. 64 -meter balloon. 
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Figure 11.- Tail-fin planforms considered for calculating tail -fin contributions to stability derivatives. 



Rolling moment of inertia, 


Density, p, kg/nr 

Figure 12.- Measured combined rolling moments of inertia of the balloon 
structure, inflation gas, and apparent mass about the axis of symmetry 
as a function of chamber density. 



Support-frame Inflation 

length, i s ,m gas 


O 9.24 Helium 

□ 9.24 Air 

O 6.24 Air 

A 3.17 Air 



+ m s<*Q- a > 2+I zz, 




Density, p, kg/m 
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Figure 13.- Measured combined yawing moments of inertia of the balloon structure, 
inflation gas, and apparent mass about the knife-edge support as a function of 
vacuum -chamber density. (Each symbol represents the average of several 
test values.) 
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